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Abstract 

Two dimensional statistical integrable models, such as the Ising model, the chiral Potts 
model and the Belavin model, becomes integrable. Because of the SU{2) symmetry of 
these models, these models become integrable. The integral models are often parameterized 
by the elliptic function or the elliptic theta function. In this paper, we study the Ising 
model, and we show that the Yang-Baxter equation of the Ising model can be written as 
the integrability condition of SU{2), and we gives the natural explanation why the Ising 
model can be parameterized by the elliptic function by connecting the spherical trigonometry 
relation and the elliptic function. Addition formula of the elliptic function is the secret of 
the exact solvability of the Ising model. 



1 Introduction 

There are many two dimensional integrable statistical models, which are classified into the 
spin model, the vertex model and the face model [1]. 

First and famous integrable and exactly solvable model is the Ising model[2]. The Ising 
model is the 2-state model, and the generalized A-state integrable model is the chiral Potts 
model[3]. The structure of the integrability condition of the Ising model is SU{2). While 
that of the chiral Potts model is the cyclic SU{2)[A, 5]. The Boltzmann weight of the Ising 
model can be parameterized by the elliptic function. 

While the 2 x 2-state integrable vertex model is the 8- vertex model[6], and the generalized 
A X A-state integrable model is the Belavin model[7]. In the Belavin model, the structure of 
the integrability condition is the cyclic SU {2)[8, 9]. The Boltzmann weight for the 8- vertex 
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model is parameterized by the elliptic function, and that for the Belavin model by the elliptic 
theta function with characteristics. 

Baxter's hard hexagon model[10] is the first integrable 2-state face model, which is ob- 
tained from the 8-vertex model by the vertex-face correspondence, and the generalized inte- 
grable model is A^^^_-^ model[ll]. The Boltzmann weight for these models are parameterized 
by the elliptic theta function with characteristics. 

Then the origin of the integrability condition of the two dimensional statistical model 
comes from the SU{2) symmetry or the variant of SU{2) symmetry. And we have the 
elliptic representation of the Boltzmann weight in many important cases. Then we expect 
the correspondence between the SU (2) symmetry and the elliptic function. In other words, 
we expect that the symmetry of the elliptic function is the SU{2) symmetry, and the origin 
of the addition formulae of the elliptic function comes from the SU (2) group structure. This 
is just the same as the origin of the trigonometric addition formulae comes from the U{1) 
group structure. For the trigonometric function, we make the connection with the circle, 
which has U{1) symmetry. For the elliptic function, we make the connection with the surface 
of the sphere, which has SU (2) symmetry. 



2 Integrability condition and exactly solvable condi- 
tion 

Let's consider the general Yang-Baxter equation of the spin model, and we explain the 
difference between the integrability condition and the exactly solvable condition. The Yang- 
Baxter equation is given by 

v{u,)U{m,u,)v{u,) = u{u^)v{u,,m)u{us). (2.1) 

This type of Yang-Baxter equation is called the integrability condition. The meaning of the 
integrability condition is that the model has the group structure so that the model has nice 
symmetry. As the Yang-Baxter equation of this type says that the product of three group 
action for two different path gives the same group action, the group structure of the model 
is expected to be SU{2) or the variation of SU{2). 

Furthermore if the Yang-Baxter equation satisfies the difference property such as 

+ U3)V{u,) = U{u,)V{u, + «3)^M, (2.2) 
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we call this type of equation as the exactly solvable condition. This is called the exactly 
solvable condition because the non-Abelian group action can be realized by the addition of 
the parameter Ui. Then the non-Abelian model is parameterized by the Abelian parameters 
Ui, which makes the model exactly solvable. The non-linear addition formula of Ui is the key 
relation that the non-Abelian integrable model becomes exactly solvable. 

The Ising model is one of the examples of the exactly solvable model. Then we pick 
up this Ising Model and explain the parameterization of the Ising model, which has the 
difference property, in a quite explicite way. The star-triangle relation of the Ising model 
can be written in the form [4, 5] 

exp(L30-^) cxi){K2a^) cxp{Lla,^) = exp{Kia^) cxp^L^a^) exp^K^a^). (2.3) 

In the previous paper [5], we give the parameterization of the Ising model by the elliptic 
function in the form 

cosh2/s:i = — sinh2i\:i = ^I^^^ (2.4) 
cn[Ui) cn[Ui) 

cosh2L* = — -, sinh2L: = ^4^^, (2.5) 

sn{K-Ui) ' sn{K-Ui) ^ ' 

(i = 1,2,3). 

Comparing Eq.(2.1) with Eq.(2.3), we have 

U{u^ = exp(Xi(7^), V{u^ = exp(L*(7^), (2.6) 

and the Yang-Baxter equation comes to have the difference property, that is, we have the 
exactly solvable condition in the form 

V{u^)U{u^ + u^)V{u^) = U{u^)V{u^ + u^)\J{u^), (2.7) 

by the parameterization of the elliptic function. We will explain why the Ising model can 
be parameterized by the elliptic function, and the Ising model comes to have the difference 
property. 



3 Differential equation for the angle and the arc in the 
spherical triangle 

The connection between the addition formula of the elliptic function and the addition formula 
of the spherical triangle has the long history, which starts from Lagrange and Legendre[12, 13] 
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and nice review article is given by Greenhill[14]. We briefly review how to parameterize 
spherical angle and spherical arc by the elliptic function, which give the parameterization of 
the Yang-Baxter relation by the elliptic function. 



3.1 First differential equation 

In this subsection, we will give the differential equation for the angle on the sphere. Various 
spherical trigonometry formulae are given by Todhunter[17]. The proof of the necessary 
formula we use here is given in the Appendix. 
The first law of cosine is given by 

cos(ai) = cos(a2) cos(a3) + cos(Ai) sin(a2) 8111(03), (3-1) 

cos(a2) = cos(a3) cos(ai) + cos(A2) sin(a3) sin(ai), (3-2) 

cos(a3) = cos(ai) cos(a2) + cos(/l3) sin(ai) sin(a2), (3.3) 



and the law of sine is given by 

sin(Ai) sin(A2) sin(A3 



sin(ai) sin(a2) sin(a3) 



= k. (3.4) 



We keep k to be constant, and we fix 03, which means that we fix through Eq.(3.4). 
Next we differentiate Eq.(3.3) and we have 

= — sin(ai) cos{a2)dai — cos(ai) sin(a2)da2 

+ cos(A3) cos(ai) sin(a2)rfai + cos(A3) sin(ai) cos{a2)da2 
— (— sin(ai) cos(a2) + cos^A^) cos(ai) sin{a2))dai 

+(— cos(ai) sin(a2) + cos(743) sin(ai) cos(a2))cia2- (3.5) 
Substituting expressions 

, . , cos(ai) — cos(a2) cos(a3) 

cos A = ' ^ 

sm(a2j sm(a3j 

cos(a2) — cos(a3) cos(ai) 



cos{A' 
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sin(a3) sin(ai) 



cos(a3) - cos(ai) cos(a2) 

cos^a) = , 3.6) 

sm(ai) sm(a2) 
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into the first term of Eq.(3.5), we have 



(First term) = (— sin(ai) 005(02) + 005(^3) cos(ai) sm{a2))dai 
— sin^(ai) cos(a2) + cos(ai) 003(03) ~ cos^(ai) cos(a2) 



smai 

— cos(a2) + cos(ai) 003(03) 



dai 



sm oi 

= — 003(^42) sin a3ciai. (3.7) 

Similarly, for the seoond term of Eq.(3.5), we have 

(Seoond term) = (— oos(ai) sin(a2) + 008(^13) sin(ai) cos{a2))da2 

= — oos(Ai) sin 03(^02. (3.8) 

Combining Eq.(3.5), Eq.(3.7), Eq.(3.8), we have 

+^^ = 0. (3.9) 



oos(Ai) 003(^2) 

Using the relation sin(Ai) = A;sin(ai), wc have oos(Ai) = ^1 — A;^ sin^(ai) where we assume 
Ai is aoute. Then Eq.(3.9) oan be written in the form 

+ = = 0. (3.10) 



We put X = 3in(ai), y = sin(a2), then we have dx = cos{ai)dai = \/l — x^dai, whioh gives 

dx 1 1 

dai = . , , = , =. Then Eq.(3.10) is written in the form of the 

v^T^ V^l-)t2 3in2(ai) ^1^^^ ^ 
differential equation of the elliptio funotion 

dx (hi 
I ^ + ^= =0. 3.11 

^{l-x''){l-W) ^(l-?/2)(l_A;V) 

Up to now, we oonsider 03 to be fixed. For the general ease, we put z = sin(a3) and obtain 
the general differential equation of the form 

dx dy dz 
I t I — Q (^3 X2) 

^{l-X^){l-k^X^) ^(l-|/2)(l_A;V) ^J{l- Z^){l-k^z'') 
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where we assume that three angels Ai, A2, A3 are all acute. If we put 

dx 



u, = r , (3.13) 

J Jil-x^){l-kH^) 



u.^[\ (3.14) 
^(1-|/2)(1-A;V) 



Us = r (3.15) 

1*1 + 1*2 + 1*3 = Const., (3.16) 

we have x = sin(ai) = su(mi, k), that is, ai = am(Mi, k). Similarly we have 02 = am(M2, k), 
as — am(u3. A;). 

3.2 Second differential equation 

In this subsection, we will give the differential equation for the arc on the sphere. We start 
from the second law of cosine 

— cos(Ai) = 005(^2) 005(^3) — cos(ai) sin(A2) sin(743), 

— cos(A2) = 005(^3) cos(Ai) — cos(a2) sin(743) sin(74i), 

— cos(743) = cos(74i) 003(^42) — cos(a3) sin(74i) sin(742), (3-17) 

and the law of sine 

sin(ai) _ sin(a2) _ sin (03) _ 1 

sin(>li) ~ sin(^2) ~ sin(A3) ~ k' K ■ ) 

The differential equation for this system is obtained by replacing a, — )■ tt — Aj(i = 1, 2, 3) 
and A; ^ 1/A; in Eq.(3.10). Then we put X = sin(y4i), Y = sin(A2), Z = sin(A3), and we 
obtain 



^(1-X2)(1-X2/A;2) ^ y'(l-r2)(i_y2/^2) 
■ = 0, (3.19) 

1- z^){i- zyk^) 
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where we assume that three length of arc oi, 02, 03 are all acute. If we put 

dX 

Ui= , (3.20) 
•f ^(1-X2)(1-XVA;2) 

U2 = r I (3.21) 

U^-r , (3.22) 
J ^(1 -Z2)(1-ZVA;2) 

C/i + C/2 + C/3 = Const., (3.23) 

we have Ai — am(C/j,l/A;) {i — 1,2,3). We also have another expression X — sin(74j) = 
sn(?7i, 1/k) = ksn{Ui/k, k). Using the law of sine, we have A;sin(ai) = sin(Aj), and sin(ai) = 
sn(-Uj, A;), sin(Aj) = ksn{Ui/k, k), we can connect Ui and Ui through sn{ui,k) — sn{Ui/k,k), 
so that we have Ui = kui, that is, 

sin(aj) = sn(iij, k) (3.24) 
sm{Ai) — sn{kui, 1/k) — ksn{ui, k). (3.25) 

We can obtain the addition formula of the elliptic function from the addition formula of 
the spherical cosine formula. Prom Eq.(3.2) and Eq.(3.17), we have 

005(02) = 003(03) cos(ai) + 003(^2) sin(a3) sin(ai), (3.26) 
— 005(^2) = 003(^3) cos(Ai) — cos(a2) sin(/l3) sin(yli), (3.27) 

and elliminating cos(A2) or cos (02), we have 

cos(«i) cos(«:i) - siii(«i) c()s(.4i) siii(«3) cos(.43) 

cos(a2) = 7 — r Tir^ 7 — ^ nr^ > (3.28) 

1 — sm(ai) sm(Ai) sm(a2) sm(>i2) 

- 003(^2) = cos(^i) co3(A3) - co3(ai) 3in(Ai) 003(03) sin(A3) 
1 — sin(ai) sin(Ai) sin(a2) sin(A2) ' ' 

Prom here, we take A2 to be obtude to find the addition formula of the elliptic function 

in the convenient form. In this case , wc have U2 = Ui + u^, and sin(aj) = sn[ui, k), sin(y4j) = 
ksm{ai), cos(aj) = cn{ui,k) {i = 1,2,3), cos(yli) = dn(tii,A;), 003(^42) — —dn{u2,k), 
cos(yl3) = dn(M3,A;). Then Eq.(3.28) and Eq.(3.29) gives 

cn{u2, k) — cn{ui + U3, k) 
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cn{ui, k)cn{u3, k) — sn(ai, A;)dn(-Ui, k)sYi{u^, k)(h\{u^, k) 
1 — A;2sn2(Mi, A;)sn2(M3, k) ' 
dn(-U2, A;) = dn(-Ui + 1x3, k) 

dn(Mi, /c)dn(M3, k) — k'^sn{ui, k)cn{ui, k)sn{u3, k)cn{u3, k) 
1 — k'^sn^{ui, A;)sn2(M3, k) 



These are convenient form of addition formula of the eUiptic function. 



4 Parameterization by the elliptic function 

The spherical trigonometry relation such as laws of the first and the second cosine and the 
law of sines are obtained from the group theoretical relation of SU (2) in the form 

exp{i74i J3.} exp{ia2 Jz} ex^{iA^Jr^} — exp{ia3 J^} exp{i(7r — A2)Jx\ exp{iai J^}, (4.1) 

where J.^, Jy, are SU{2) generators with [J^, Jy] = iJ^, [Jy, J^] = iJ^, [J^, Jx] = iJy For 
J — 1 case, Eq.(4.1) becomes 

{1 - (1 -cos(/li)) Jl + ism{Ai)Jx}{l - (1 - cos(a2)) J| + i sin(a2) Jj 
{1 - (1 - cos(A3)) Jl + i sin(A3) = {1 - (1 - cos(a3)) + i sin(a3) 
{1 - (1 -cos(7r - A2)) + isin(7r - A2)Jx}{l - (1 - cos(ai)) + i sin(ai) J;,}, 



and using sn(ii3, k) — Jl — cn2(w3, /c), we have 



sn(M2, k) = sn(Mi + M3, k) 

sn(Mi, /c)cn(M3, k)dn{u3, k) + cn(-Ui, k)dii{ui, k)sn{u3, k) 



(3.32) 



1 — fc2sn2(-Ui, A;)sn2(-U3, k) 
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^2 

Fig.l The change of vectors on the spherical triangle 



This relation is understand as the integrabihty condition on the sphere. We consider 3 
points on the sphere Pi, P2, P3. Starting from the vector ni at point P3. Then the vector 
parallel transports along the tangential direction of the large circle or rotates around the 
point Pi{i = 1,2,3). The following group element corresponds to the parallel transport or 
the rotation around the point Pj(i = 1, 2, 3) of the vector in the following way 

exp{iA3Jx} : ni 115, ex.-p{ia2Jz} ■ ns ne, exp{iAi J^,} : ne 114, 
exp{iai J2} : ni n2, exp{i(7r - A2)Jx} ■ ^2 na, expfiaaJ^} : ns ^ n4 

Therefore Eq.(4.1) is the integrabihty condition that two different ways to change the vector 
from ni to n^, that is, path l:{ni — 7> ns — ne ^ n^} and path 2:{ni — )■ n2 — )■ na — )■ n4}, 
gives the same group action, which is nothing but the integrabihty condition of SU{2). 
While the star-triangle relation can be written in the form 



exp{2L* J^} exp{2K2Jz} exp{2L* J^} = exp{2is:3 Jj exp{2L; J^} exp{2is:i Jj. (4.2) 



In the previous paper, we parameterize the Ising model with the elliptic function in the form 



cosh{2Ki) 



cosh(2L*) 



cn(vi, k') ' 
dn(vi, k') 



1 



sinh(2L*) 



sn(f;, k') 
cn(t'j, k') ' 
/csn(vi, k') 



(4.3) 



(4.4) 



cn(vj, k') ' 



cn{vi,k') ' 
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V2^vi + v^, (i = l,2,3). (4.5) 
By using the Jacobi's imaginary transformation, we put Vi = iui and we have 

sn{v,,k')^^'-^pl^^, cn(^„fcO = ^^, dn(^„ fc') = (4-6) 
cn(tij, k) cn[Ui, k) cn[Ui, k) 



then we have 



Using the relation 



cosh(2i<'j) = cn(ui,k), sinh(2i<'i) = isn{ui,k), (4.7) 
cosh(2L*) = dn{ui, k), sinh(2L*) = iksn{ui, k), (4.8) 
U2^ui + us, (i = l,2,3). (4.9) 



sn{u, 1/k) — ksn{u/k, k), cn{u, 1/k) — dn{u/k, k), 

dn{u, 1/k) = cn{u/k, k), (4.10) 



and putting Ki = iKi, L* = iL*, we have 



cos(2i^i) = cn(?7,,, A;), sin(2i^i) = sn{ui, k), (4.11) 
cos(2L*) = cn(A;Mi, 1/k), sin(2L*) = sn{kui, 1/k), (4.12) 
U2^ui + U3, (i = 1,2,3). (4.13) 



For the spherical triangle, we take A2 as obtuse, and Ai and A3 are acute. In this case 

dA2 dAi dA^ 



l-sm\A2)/k^ Jl-sm\A,)/k^ Jl-sm\A3)/k 



(4.14) 



and we take U2 = Ui + U3, which gives U2 — Ui + U3 by using C/j = kui{i — 1, 2, 3). While 
the spherical triangle is parameterized by 

cos(ai) = cn(itj, /c), sm{ai) — sn{ui,k), (4-15) 

cos(Aj) = cn{Ui, 1/k) — cn{kui, 1/k), sin(Aj) = sn{kui, 1/k), (4-16) 

U2^ui + U3, {i^ 1,2,3). (4.17) 

Then we have 

2Ki = = a,m{ui, k), 2L* = A = a,m{kui, 1/k), (4-18) 
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where we must notice that A2 is obtuse, so that tt — ^2 is acute. Then the Yang-Baxter 
equation is written in the form 

exjp{iam{kui, l/k)Jx} exp{ieim{u2,k)Jz} ex.p{iam{kus, l/k)Jx} 

— exp{iam{u3, k)Jz} exp{isim{ku2, J^} exp{mm(ui, A;) J^}. (4.19) 



5 Abel's addition theorem 

The addition theorem of the algebraic function, that is, the Abel's addition theorem is the 
essential key relation to solve the statistical model. Then we consider the old problem of 
the Abel's addition theorem here. We put /^{x) = Aq + Aix + ^422;^ + ^30;^ + A4X^, and we 
consider the g — 1 algebraic curve = f4{x): then Jacobi or Richelot[15, 16] tried to find 
the conserved quantity of the differential equation 

+^£^ = 0. (5.1) 



While we consider here the Abel's addition theorem of the standard type, that is, 

^"^^ +^£^ + ^£^ = 0. (5.2) 



^Jh{xi) \IUx2) ^Jh{xz) 

In general, we take 

f2n-2{x) =Ao + AiX + A2X^ + ■ ■ ■ + A2n-iX^'^~^ + A2n-2X^"'~^ , (5.3) 

and consider the differential equation of the form 



dxi dx2 dxn 

+ , = + --- + ^ =0, 



\lf2n-2{xi) ^Jhn-2{X2) hn-2{Xn) 



XidXi ^ X2dX2 ^ ^ XndXn _ ^ 



\lf2n-2{xi) \lf2n-2{x2) f2n-2{Xn) 

Xi^'^dXi X2'^'^dX2 , , Xn"-^dXn „ .s 

y + , H \- I = 0. (5.4) 

\lf2n-2{xi) \]f2n-2{x2) f2n-2{Xn) 
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We define F{x) by 

F{x) = {X - Xi){x -X2)---{X - Xn-l){x - Xn), (5.5) 

then we have 

F'{xi) = {xi - X2){xi - X3) • • • {xi - Xn-i)ixi - (5.6) 
etc. Then if Xi, which depend on t, satisfies the following equation 

dXi _ Xi^J f2n-2{Xi) 



dt F'{xi) ' 

we can show that Eq.(5.4) is satisfied. Here we use the following theorem for F{x) 



(5.7) 



_^ = V— 7^-^^, (A; = 0,1, •••,71-1). (5.8) 

F{x) {^^F'{xi){x-Xiy ^ ' ' ' ^ y ) 

The proof is the foUowings: As F{x) is the n-th polynomial and it is much higher polynomial 
than x^ , we can write 

x^ " a. 



^E? —V (A) = 0,l,^^^,n-1), (5.9) 



then we multiply {x — Xj) and take the limit x ^ Xj, we have 
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= E'^»^jj = %' (5.10) 



which gives Eq.(5.8). By multiplying x and taking a; — >■ oo in Eq.(5.8), we have 

t.-E^^=^k,n-i, ik^O,l,---,n-l). (5.11) 

i=l ^ K-'^i) 

Then the differential equation of the problem is satisfied in the following way 

^ X.'^dx, ^ ^ X,'^+'^f2n-2{x^)dt ^ ^ X.^+^rft ^ ^ 
i=l \Jf2n-2{Xi) i=l F'{Xi)^ f2n-2{Xi) i=l F' (Xi) 

{k = 0,l,---,n-3) 
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5.1 Abel's addition theorem I: one conserved quantity 

We will integrate this differential equation and find the conserved quantity. Starting from 
Eq.(5.4), we have 



d^Xi _ dXi d ( Xi^J f2n-2{,Xi)\ dX2 d ( Xl^J f2n-2{xi) 

~d^ " ~~di~dxiy F'{xi) J ^ "dTda^ \^ F'{xi) 

dXn d ( Xl^f2n-2{xiy 

^~dt~dx'n I F'{xi) 



+ 



(5.13) 



XiyJ f2n-2{xi) d ( XisJ f2n-2{xi)\ XiX2^J f2n-2{xi) f2n-2{x2) 1 

F'(xi) I F'(xi) ) ^ F'{X^)F'{X2) (Xi - X2) ^ 



y f2n-2{Xl) f2n-2{Xn) 1 , 

+ F'{x,)F'{xn) (xi - x„) ^^-^ ^ 

^1 d ( xlf2n~2{Xi) \ ^ XiXi^f2n-2{Xl)f2n-2{Xi) 1 

2(ixi 1, (F'(xi))2 F'(xi)F'(x,) {x.-XiY ^ ■ ^ 

By taking the sum, we have 

While we use the formula 

X^f2n-2ix) -\^( ^'ihri-2{Xi) d f X^ f2n-2{Xi) \ 1 ^/c^17^ 

F{xf ^ \F'{x,f{x - x,y ^ dx, \ F'{x,Y I {x - X,) ) 

The proof of the above is the following: We can write 

^^hn-2{x) aj hi \ 

F{xY - - h ■ ^ ^ ^ 

We multiply {x — XjY and take the limit x — > x^-, the we have 

x']f2n-2{.Xj) 



F'{x,y 



^aAj = %. (5.19) 
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Next we multiply {x — xj)^ and further multiply — and take the limit x ^ xj, we have 

(Jjtjb 



dx F x 2 dxi F' x,y h ^ 



d {X- Xjfx^f2n-2{x) _ d x]f2n-2 {Xj) _ ^ 

Cj F'{xj 

From Eq.(5.18), Eq.(5.19), Eq.(5.20), we have Eq.(5.17). While by using the relation 

X^f2n-2{X 



Fix) 
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A2n-2 = [A^n-^ + 2^2^-2 jZ^^j^ + O , (5.21) 



and by comparing the coefficient of 0{—) in Eq.(5.17), we have 

+ 2A2n-2 E a;. = V f a;f /a^.^ (x,) \ ^^^j^ ^ ^^^^ 

tidxi\ F'{xiY J dt^h 

n 

where we use Eq.(5.16). We put p = ^^Xi, then we have the differential equation of the 

i=l 

form 

2^ = A2n-3 + 2A2n-2P. (5.23) 

dp 

Multiplying — in both side of Eq.(5.23), we have 
dt 

Then we have one conserved quantity 

E ^^V^l^) = A2^.s tx. + A2^^2 itxX + C,. (5.25) 



1=1 \i=l 



5.2 Abel's addition theorem II: another conserved quantity 

As Richelot[16] showed, wc transform the original differential equation of Xi into that of 
= — . In this variable, we have 

Xi 



dxi = (5.26) 

Si 
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1 1 

, = ; = , (5.27) 

^f2n-2{Xi) yjhn-2{l/^i) V^2n-2(Cj) 

where 



92n-2{^i) = ^2n-2 + A2„_3e + ' ' " + ^iC'""' + ^0^'""', (5.28) 



which gives 



^ ' (A; = 0,l,---,n-3) 



\J hn-2{Xi) 



(A;' = n-3-A; = 0,l,---,n-3). (5.29) 



\Jg2n-2{^i) 

Therefore the differential equation is transformed into the form 



Y: il^ = 0, (A; = 0,l,---,n-3). (5.30) 

Yfl'2n-2(6) 



i=l 



Defining Fi(^) = nr=i(C ~ C?) ^^'^ comparing with Eq.(5.25), we have another conserved 
quantity 



Using 



E^S^i =AY.(.+Mj:(,r+c,. (6.31) 

=1 ^li^i) / i=l i=l 



m^) = ^ (5.32) 



we have 



E (^1^2 • • • :r„)^ -A±^ + M± + C2. (5.33) 



=1 -^i^ I i=l i=l 
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5.3 Simple n = 3 case: Abel's addition theorem for the elliptic 
function 



For n — 3 case with f4{x) = (1 — — k'^x^), we numerically checked by Maxima the 

following addition formula 



6 Summary and discussion 

The integrability condition, which is called the Yang-Baxter equation, in two dimensional 
statistical models means that such integrable models have the group structure. The Boltz- 
mann weight of such integrable models often can be parameterized by the elliptic function 
or the elliptic theta function. Furthermore if the Yang-Baxter equation have the nice dif- 
ference property, we can exactly solve the model. Such exactly solvable models often can 
be parameterized by the eUiptic function. If we notice that the Yang-Baxter relation is the 
relation of the products of the three group action, we expect that the integrable model has 
the SU(2) or SU{2) variant group structure. Furthermore if the model is exactly solvable, 
we expect that the exactly solvable model can be parameterized by the elliptic function. As 
the Ising model is the exactly solvable model, the above expectation is realized by rewriting 
the Yang-Baxter relation in the form of the integrability condition of SU{2) group. Further- 
more we can parameterize the Boltzmann weight by the elliptic function, which comes to 
have the difference property, by making the connection between the SU{2) group element 
and the elliptic function. 

In this way, we can understand the reason why the Ising model can be parameterized by 
the elliptic function in a quite natural way, where the connection between the SU{2) and 
the eUiptic function is essential. Addition formula of the eUiptic function is the secret of the 
exact solvability of the Ising model. 

Finally, we find some of the conserved quantity for the hyperelliptic differential equation, 
which gives the explicit form of the Abel's addition theorem. This will will be useful to find 




(5.34) 



(5.35) 



where 

Xi = sn(tij). 
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the group structure of such hypereUiptic function. 
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A Spherical trigonometry formulae 



In this appendix, we give the brief proof of various spherical trigonometry relations. We 
consider the sphere with unit radius and we put three points Pi, P2, P3 on the sphere. We 
consider three unit vectors which go from the origin to points Pi, P2, P3, and we denote such 
unit vecotors as rii, 112, ns . Further, we denote the length of the arc connecting and 
as Oi, and that connecting ng and iii as 02, and that connecting rii and 112 as 03 . 




Fig. 2 Angles and circs of the spherical triangle 



Then we have 



cos(ai) = n2 • na, cos(a2) = na • ni, cos(a3) = ni • n2. 



(A.l) 



We make the dual unit vectors perpendicular to ni and nj = 1,2,3), {i 7^ j) in the 
following way 



n 



n2 X ns 



|n2 X nsl |n3 x ni| |ni x n2| 



ns X ni 



ni X n2 



(A.2) 



For planes perpendicular to ni and nj, we denote = 1,2,3), {i ^ j), which are given 

by 



^23 : Hi • X = 0, Zsi : n2 • X = 0, Zi2 : ng • x = 0. 



(A.3) 
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We denote the rotation angle of the spherical triangle around vector rii as Ai as is shown in 
Fig.2, then we have cos(Ai — tt) = iij • 113 from geometrical consideration. Similarly we have 

cos(Ai — n) — 112 • 113, cos(742 — tt) = 113 • n^, cos{A3 — n) — nl • 112. (A.4) 

Using the formula 

(a X b) X (c X d) = |a, c, d|b — |b, c, d|a = |a, b, d|c — |a, b, c|d, (A. 5) 

where 

|a, b, c| = a • (b X c) = (a X b) • c = (determinant of a, b, c), 

we have 

ni = ^ fj, n2 = ^ ^, n3 = ^ ^. A.6 

X njl |n5 x ni\ \nl x 

A.l Law of cosines 

Using the theorem 

(a X b) • (c X d) = (a • c)(b • d) - (a • d)(b • c), (A.7) 

we have 

/ 4 ^ * * (ns X m) (ni x na) 

cosfTT — Ai) — nX nl — ■ f 

|n3 X nil |ni x na] 

(rii ■ n3)(ni ■ na) - (n2 ■ n^) cos(a2) cos(a3) - cos(ai) 



|n3 X ni||ni x na] sin(a2) sin(a3) 

Then we have the first law of cosine in the form 



(A.8) 



cos(ai) = cos(a2) cos(a3) + cos(Ai) sin(a2) sin(a3), 
cos(a2) = cos(a3) cos(ai) + cos(A2) sin(a3) sin(ai), 
cos(a3) = cos(ai) cos(a2) + cos(743) sin(ai) sin(a2). (A.9) 

While from the following relation 

(n* X n*) (n^ x n*) _ (n^ • n*)(n* • n*) - (n* • n*) 



cos(ai) = n2 • na 



\nl X nJI |nj x rial \nl x nJUni x rial 



cos(7r - A2) cos(7r - A3) - cos(7r - Ai) 
sin(7r — A2) sin(7r — ^13) 
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Then we have the second law of cosines in the following form 



— cos(a4i) = cos(y42) cos{As) — cos(ai) sin(A2) sin(A3), 

— cos(A2) = 005(^3) cos(Ai) — 005(02) sin(A3) sin(Ai), 

— 005(^3) = cos(Ai) cos(y42) — 005(03) sin(y4i) sin(A2). (^-11) 

A. 2 Law of sines 

Using the relation 



(na X nij (ni x naj 





n3,ni,n2| 




na X ni 


|ni X 



n3 X rill |ni x 

ni = rill, (A.12) 





ni,n2,n3 




na X ni 


ni X n2| 



we have 



X Hgl \nl X n^l x 



n2 X na na x ni rii x n2 





ni,n2,na 




|n^n*,n*| 


|ni X n2| 


n2 X na 


na X ni 




ni,n2,na| 



(A.13) 



where we use 

1***1 /* *\* -^1 ^ -^2 •■ -^3 / \ 

Hj, Hg = (n^ X Ha) • Hg = -ng • (ni x n2) 

|n2 X nallna X ni||ni X n2| 

(|ni.n2,na|)^ 

|n2 X nallng x ni||ni x n2| 

From Eq.(A.13), we have the law of sines in the following form 

sin(7r — Ai) sin(7r — A2) sin(7r — As) 



(A.14) 



sin(ai) sin(a2) sin(a3) 



or 



sin(Ai) sin(A2) sin(A3) 
sin(ai) sin(a2) sin(a3) 

We use Eq.(A.9), Eq.(A.ll), Eq.(A.15) in the body of the text. 



(A.15) 
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